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Confidence intervals and joint confidence sets are constructed for 
tlie nonparametric calibration of exponential Levy models based on 
prices of European options. This is done by showing joint asymptotic 
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CN 1 Introduction 

O 

The unknown future development of financial markets faced by its partici- 
pants, whether as investors, as traders or as companies, can be understood 
•'H to consist of model risk and "Knightian uncertainty" |10l I17j . The first 

describes the risk for a given calibrated model and can be evaluated by 
probabilistic methods, whereas the second incorporates the lack of knowl- 
edge on the underlying probability measure and is typically treated by worst 
case scenarios, for example, by stress testing, which amounts to taking the 
supremum or infimum over a range of probability measures. 
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By the choice of the model, there is a trade-off in the uncertainty between 
the cahbration error and the misspecification of the model, where larger 
models reduce the latter. Since the calibration error is statistically traceable 
and the misspecification is hard to track, we choose a rich model such as 
a nonparametric one and then asses the calibration error by constructing 
confidence sets. 

More precisely, we consider the nonparametric calibration when the risk- 
neutral price of a stock {St) follows an exponential Levy model 

St = 5e''*+^* with a Levy process for t > 0. (1.1) 

A thorough discussion of this model is given in the monograph by Cont and 
Tankov [7]. They introduced in [8l|9] a nonparametric calibration method for 
this model based on prices of European call and put options, in which a least 
squares approach is penalized by relative entropy. Belomestny and Reifi [2] 
used a different approach to the same estimation problem, where they reg- 
ularized by a spectral cut-off and constructed estimators that achieve the 
minimax rates of convergence. We show asymptotic normality and construct 
confidence sets and intervals for their estimation procedure. Methods simi- 
lar to theirs were also applied by Belomestny p,j to estimate the fractional 
order of regular Levy processes of exponential type and by Trabs [22j to 
estimate self-decomposable Levy processes. 

Confidence sets measure how reliable the estimation is. This is particu- 
larly important if the calibrated model is to be used for pricing and hedg- 
ing. For a recent review on pricing and hedging in exponential Levy models 
see |21] and the references therein. For the influence of model uncertainty 
on the pricing see [6]. 

Nonparametric confidence intervals and sets for Levy triplets have not 
been studied with the notable exception of the work by Figueroa-Lopez |12] . 
The work is more general in the sense that beyond pointwise confidence inter- 
vals also confidence bands are constructed. On the other hand, the method 
is based on high-frequency observations so that the statistical problem of 
estimating the Levy density reduces to the problem of density estimation. 
In contrast to this direct observation scheme our method is based on option 
prices and thus the calibration is a nonlinear inverse problem, which is mildly 
ill-posed for volatility zero and severely ill-posed for positive volatility. 

Confidence intervals and sets in nonparametric problems are a subtle is- 
sue. Whether adaptive confidence intervals for the estimators of the volatil- 
ity, the drift and the intensity exist is an open question. We show asymptotic 
normality for the parametric estimators of the volatility, the drift and the 
intensity. We also proof asymptotic normality for the pointwise estimators 
of the Levy density. The joint asymptotic distribution of these estimators 
is derived in both the mildly and the severely ill-posed case. This is used 
for the construction of confidence intervals and joint confidence sets as well 
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as for a test on the value of the volatihty. The asymptotic normahty results 
are based on undersmoothing and on a linearization of the stochastic errors. 

The paper is organized as follows. The model and the estimation method 
are given in Section [2] The main results are formulated in Section |3j They 
are applied to confidence intervals and to a hypotheses test on the value of 
the volatility in Section |4j We conclude in Section [5} Proofs are deferred to 
Section [6] and Section [7| Uniform convergence is treated in Section [8| 



2 The model and the estimators 

We denote by C{K, T) and V{K, T) prices of European call and put options 
on the underlying (St) with strike price K and maturity T. We suppose 
that the risk-neutral price of the stock (St) follows the exponential Levy 



model (1.1) with respect to an equivalent martingale measure Q, where 
5 > is the present value of the stock and r > is the riskless interest rate. 
We fix some T and assume that the observations are given by option prices 
for different maturities (Kj) corrupted by noise: 

Yj =C{Kj,T) + ajej, j = l,...,n. 

The noise levels (cjj) are assumed to be positive and known. The minimax 
result in [2] is shown for general errors (e^) which are independent centered 
random variables with E[e^] = 1 and supjE[ej] < oo. We transform the 
observations to a regression problem on the function 



0{x) :-- 



S-^C{x,T), x>0, 
S-^V{x,T), x<0, 



where x := log{K/S) — rT denotes the log- forward moneyness. The regres- 
sion model may then be written as 

Oj = 0{xj) + 5j€j. (2.1) 

We call the volatility of a Levy process a, the drift 7 and the inten- 
sity A. As in |2j we consider only Levy processes (Xt) with a jump com- 
ponent of finite intensity and absolutely continuous jump distribution. For 
a jump density v we denote by iJ,{x) := e^i^(x) the corresponding expo- 
nentially weighted jump density. The aim is to estimate the Levy triplet 
T = (a"^,7,/i). In the remainder of this section we present the spectral 
calibration method of Belomestny and Reii^ [2J . The method is based on a 
pricing formula by Carr and Madan [5] , which relates the Fourier transform 
FO{u) := 0(x)e*"^dx to the characteristic function ipriu) := E[e*""'^^]. 
That is why we define 



tl^iu) := ^ log (1 + iu{l + iu)FO{u)) = ^ \og{^pT{u - i)) 



(2.2) 



+ i(o-2 + 7)ti + (0-V2 -h 7 - A) + -F/x(n), 
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where the first equahty is given by the above mentioned pricing formula and 
the second by the Levy-Khintchine representation. This equation links the 
observation of O to the Levy triplet that we want to estimate. Let be 
an approximation on the true function O. For example, can be obtained 
by linear interpolation of the data (2.1). We further define the empirical 
counterpart of ^p by 

Aiu) ■■= ^ log>«,(„) (1 + iu{l + iu)J='0,{u)) , 
where the trimmed logarithm log>^ : C\{0} — )• C is given by 

W (z) - / - 

The logarithms are taken in such a way that tp and are continuous with 
V'(O) = il^eiO) = and k{u) G (0,1) is specified in [2j. Considering ( [2^ 
as a quadratic polynomial in u disturbed by J-'fj, motivates the following 
definition of the estimators for a cut-off value U > 0: 



a' 



u 



rU 

7 := -CT^ + / lm{ilj^{u))w!:! {u)du, 
J-U 



A := 



+ 7- 



u 



(2.3) 
(2.4) 
(2.5) 



where the weight functions w^^ , and satisfy 



u 



-wl^{u)du = 1, 



uw^ {u)du = 1, 



u 



Wxiu)du = 1, 



(n)du = 0, 



u 



u 



u'^w'^{u)du = 0. 



(2.6) 



The estimator for /i is defined by a smoothed inverse Fourier transform of 
the remainder 



fi{x) := J" 



-1 



'^e{u) + —{u-i) - i^{u- i) + X] w^{u) 



(x). (2.7) 



The weight functions for all [/ > can be obtained from w^, w}^, w\ and 
wj^ by rescaling: 

w^iu) = U-^wliu/U), wl^iu) = U-^wl^{u/U), 
w'iiu) = U-W,{u/U), w^i (n) = wliu/U). 
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Since ipe{—u) = V'e(^) only the symmetric part of w^, w\ and the anti- 
symmetric part of w}^ matter. The antisymmetric part of wj^ contributes a 
purely imaginary part to Without loss of generality we will always 

assume w^, w\, wj^ to be symmetric and u;^ to be antisymmetric. We fur- 
ther assume that the support of w^, w}:^, w\ and wj^ is contained in [—1,1]. 
We define the estimation error and likewise for the other 

estimators. We will also use the notation Aip^ := — ijj. The estimation 
error Aa^ can be decomposed as 



Act^ : = 2[7-2 f Yle{Tfi{Uu))wl{u)du 
Jo 

Jo 



{21 



The first term is the approximation error and decreases in U due to the decay 
of J-^. The second is the stochastic error and increases in U by the growth 
of A^/^g. The cut-off value U is the crucial tuning parameter in this method 
and allows a trade-off between the error terms. The other estimation errors 
allow similar decompositions as Ao"^ in (2.8). 



We shall analyze the asymptotic properties of the stochastic errors in 
depth. To bound the approximation errors some smoothness assumption is 
necessary. We assume that the Levy triplet belongs to a smoothness class 
Qs{R-, Cmax) with s G N and i?, fimax > specified in [2t Definition 4.1]. The 
assumption T G Gs{R-,(^raa.yi) includes a smoothness assumption of order s 
on /i leading to a decay of F^. To profit from this decay when bounding 
the approximation error, we assume the weight functions to be of order s, 
this means 

T{wl{u)/u'),F{wl^{u)/u'),F{wl^{u)/u'),F ((1 - wl{u))/u') G L\n). 



(2.9) 

3 Asymptotic normality 
3.1 The main results 

The aim of this section is to establish asymptotic normality results for the 
estimators. We would like to state that the appropriately scaled errors of 
the estimators converge to normal random variables. The starting point 



of the error analysis is the decomposition (2.8) into the approximation er- 
ror and the stochastic error. The approximation error is deterministic and 
only the stochastic error can be expected to converge with appropriate scal- 
ing to a normal random variable. It is common practice to resolve this 
problem by undersmoothing, which means that the tuning parameter is 
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chosen such that the approximation error becomes asymptotically negli- 
gible. Thus the cut-off value has to grow fast enough. This is ensured 
by the condition eC/(e)^^*"^^)/^ — t- oo in the case of volatility zero and by 
e[/(e)''"'"^ exp(Ta"^C/(e)^/2) — )• oo in the case of positive volatility. 

Since the approximation errors are negligible by these conditions we 
focus on the stochastic errors. To simplify the asymptotic analysis we do 



not work with the regression model (2.1 ) but with the Gaussian white noise 



model. This is an idealized observation scheme, where the terms are easier 
to analyze. At the same time asymptotic results may be transferred to the 
regression model. The Gaussian white noise model is given by 

dO,{x) = 0{x)dx + e6{x)dW{x), (3.1) 

where is a two-sided Brownian motion, 6 £ L^(]R) and e > 0. Here 
the empirical counterpart of O is given directly by the model and no 
further approximation is necessary. Transferring asymptotic results from 
the Gaussian white noise model to the regression model is formally justified 
by the concept of asymptotic equivalence. Due to the asymptotic equivalence 
between regression with regular errors and regression with Gaussian errors 
in |14j we will consider for our asymptotic analysis only the case of Gaussian 



errors. The asymptotic equivalence of regression (2.1) with Gaussian errors 



and the Gaussian white noise model (3.1) is given in Brown and Low [4 



where 6j = 5{xj) and e corresponds to l/\/n up to a logarithmic factor. 
Further assumptions for the asymptotic equivalence to hold are given in 
Section [6l 

The stochastic errors involve the term A'iIj^{Uu), which is a difference 
between two logarithms. For z,z' £ C\{0} and k > it holds log>f.(2) — 
log(z') = log>^/|^,| {z/z'). That yields 



oo 

iUux i 



where {u) := k{Uu)/\1 + iUu{l + iUu)J'0{Uu)\ < 1/2, see ^ (6.3)]. We 
define a linearization C^^u of the logarithm and the remainder term TZe^jj by 

^-'C' - T(i..^:';^;^(^„)) £ ^""m^i (3.2) 

7^e,{/(u) := Aij,{Uu) - C,^u{u). (3.3) 

To ensure continuity of the Gaussian process X{u) = e^^^5{x)dW{x) 
we assume that there is ap > 1 such that f^^{l + \x\yS{x)'^dx < oo. In [19] 
it is shown that on this assumption X satisfies the Kolmogorov-Chentsov 
criterion [.16l p. 57] and thus has a continuous version. In the sequel we are 
always working with this version. 
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The remainder term TZ^u in (3.3) is small when the argument of the 



logarithm is close to one, that is when Ce^u is small. Since we are integrating 



over the unit interval in (2.8) we want C^^^u to be uniformly small. We shall 
use the notation A{x) < B{x) as x — )• oo synonymously with the Landau 
notation A{x) = 0{B{x)) as x ^ oo, meaning that there exist M > and 
xo G H such that A{x) < MB{x) for all x > xq. 

Proposition 1. For all q > 1 holds 

1/9 



E 



sup \£.e,u 
«e[-i,i] 



(n)l^ 



<eC/Vlog(C/) exp{Ta^U^/2) 



as U 



oo. 



This proposition is proved in Section [6] by metric entropy arguments. In 
the following theorems we control the supremum oi C^^u and thus the remain- 
der term TZ^^u by the condition eU [e^ y^log{U (e)) exp{Ta'^U{e)^ /2) 0. 
Then the asymptotic distribution of the stochastic errors Ail>^{U u)w{u)du 
is governed by the linearized stochastic errors £e,u{u)w{u)du and the re- 
mainder term TZe^uiu)w{u)du is asymptotically negligible. In the case 
(7 = the stronger condition e[/(e)^/^ — )• is assumed, which is needed for 
the stochastic errors to converge to zero. 

In the results on asymptotic normality we will also include the estimator 
/i(0) of the jump density at zero. This only makes sense by our smoothness 
assumption on n since there is no way of detecting jumps of height zero. 
Unlike for points x 7^ it will turn out that not the weight function wj^ 
determines the asymptotic distribution but the effective weight function 

wo{u) := wj^{u) + wl{u) / v'^wj^{v)dv/2 — w\{u) / wj^{v)dv. 



J-i J-i 

The first theorem states the joint asymptotic normality result for the 
mildly ill-posed case of volatility zero. 

Theorem 1. Let a = 0. Let 6 be continuous at T7, xi + T7, . . . ,Xn + T'j 
and let Tb'^ G L^{^. For j = l,...,n let xj S 1R\{0} be distinct and 
let Vb, Wo, Wxi . . . , Wx„ be independent Brownian motions. If ell (e)^/^ — )■ 



and e?7(e)(2«+5)/2 

/ ?7(e)+i/2 
?7(e)-i/2 
C/(e)-3/2 
C/(e)-5/2 
C/(e)-5/2 

N-5/2 



— )• 00 as e 

A7 
AA 

AA(0) 
A/i(xi) 



0, then it holds 

( d(0) jlu^wl{u)dWo{u) 
d(0) u^w}^{u)dVo{u) 
d(0) u'w\{u)dWo{u) 

-> d(0) j^u'wo{u)dWo{u)l{2Ti) 

d{xi) /o\2^i(n)dI^,,(n)/(27r) 



V ^(e) 5/2 A/i(x„) y \d{xn) Jiu^wliu)dWxM/{'^^) J 
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as e — )• 0, where d{x) := 2^/^^5{x + T7) exp(r(A — 7))/T. 



Next we consider the case ct > 0. Let 6 be in L°°(]R) and ll^llLa^jj^) > 0. 
We set 

d := V2\\5h^u) exp(r(A - 7 - a^/2))T-'a-\ (3.4) 

and define real- valued random _variables VF^ u and u by u + (7 := 



2d ^ £e,i/('u)du. By Lemma 3 below 



1 / W,,u \ d f W 



eexp(Ta2f/2/2) V 7 V ^ 



(3.5) 



as ?7 — )• 00, where W and F are independent standard normal random 
variables. 

The following theorem treats the stochastic errors in the case of positive 
volatility. Since the theorem contains no statement on the approximation 



errors, the condition (2.9) on the order of the weight functions may be 
omitted. 

Theorem 2. Let a > and 5 £ L°°(1R). Assume for the cut-off value 
U{e) 00 and eU (e)^ ^/log{U {e)) e^p{Ta^U{e)^ /2) ^ as e ^ 0. Let 
^ : [0, 1] — )• IR be Riemann-integrable, in L°^([0, 1]) and contin- 



nous at one. For x it holds with the notation U := U{e) 

( 2j^Re{/\ij,{Uu))wl{u)du-dwl{l)W,^u \ 
2 /(J Im(AV'e(f/n))u'^(n)dn - dw]{l) V,,u 
2 Re{^^^{Uu))w\{u)du - dwl{l)W,^u 
\ T-^ {^i,,{Vn)wl(u)\ (Ux) - dwlil) zIu{x)/{27t) J 



0, 



as e — )• 0, where Z^^u{x) := cos{Ux)We,u + s\\i{Ux)Ve^u- 

The assumption T £ Gs{R,(^ma.x) includes a G [0,(Tmax]- The condi- 
tion e;7(e) Vlog(f^(e)) exp{Ta^U{ef /2) ^ is especially fulfilled if U{e) < 
a'^-'^(21og(e^"'^)/T)^/^ for any fj > CTmax- For the estimation it suffices to know 
some upper bound Umax of a. The theorem shows that regardless whether 
one undersmoothes or not the stochastic errors converge with appropriate 
scaling to normal random variables. For the statement on asymptotic nor- 
mality we have to undersmooth and further knowledge on the volatility is 
necessary. 

In many situations the volatility a is known or can be estimated eas- 
ily. It can, for example, be determined from high frequency data of the 
underlying asset since the volatility is the same for equivalent measures. 
In the following we will assume either that the volatility a is known as 
in Cont and Tankov [8] or that we have a sufficiently good estimator of 
the volatility. To control the remainder term we choose U{e) such that 
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e?7(e) Vlog(f/(e)) exp{Ta'^U{ef /2) ^ as e ^ 0. We also assume the 
undersmoothing condition eC/(e)*^"^ exp(Tc7^C/(e)^/2) — )■ oo as e — 0. A 
smoothness parameter s > 2 is imphcitly assumed so that both condition 
can be satisfied simultaneously. A possible choice of U (e) is 



U{e) :- 



log 



(3.6) 



where a G (1, (s + l)/2). Then it holds 

{oo f3 > 2a 
(^)" /3 = 2a 
(3 <2a 

as e — )• 0. Especially the term diverges for /3 = s + 1 and converges to zero 
for /3 G (2, 2a) so that both conditions on U{e) are fulfilled. 

Next we state the joint asymptotic normality result for the severely ill- 
posed case of positive volatility. 

Theorem 3. Let a > and let the cut-off value U{e) he chosen such that 
eU{ef^/kfg{U{e))exv{T(T^U{ef/2) and eU{eY+^ ex.-p{Ta'^U{ef /2) 
oo as e^O. Let 6 £ L°°(1R). It holds 



geT<x2[/{e)2/2 



/ U{e) 
U{e) 

U{e) 
\ U{e) 



A-y 

AX 
AA(0) 
Af{x) 



dw}^{l)V^^u{e) 

d^(l)H^.,C/{e) 

dwo{l)W,^U(e)/{2n) 

dwj,il)Z,^Uie)ix)/{27T) J 



0, 



as e — )• 0, where x G 1R\{0}, Z^^u[x) := cos(C/a;)W£^(/ -|-sin([/x)V'e,(7 and d is 



given by (3.4). 



3.2 Discussion of the results 

Theorems [l] and [s] include the asymptotic distribution of o"^ , which may be 
used for testing the hypotheses Hq : a = ao. If o" is known, we can set 
Then the statements of the theorems hold with w]^ constant to 
zero. The estimation method can give negative values for a^, A and i>(x). 
By a postprocessing step the estimated values can be corrected to be non- 
negative. 

In Theorem [T] the noise level 5 enters only locally into the asymptotic 
variance whereas in Theorems [2] and [3] the asymptotic variance depends on 
the L^-norm of 5 through the factor d. In fact for o" = it is possible to 
estimate 7 and A directly from local properties of the option function O at 
7r as remarked in [2J. This local dependence on the noise level resembles 
some similarity to deconvolution, for instance, to the case of ordinary smooth 
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error density, where the density of the observations enters locally into the 
asymptotic variance |1H [23] . For the weight functions the local and global 
dependence is vice versa. In Theorem [l] the weight functions w^, w}^, wl, 
Wq and wji^ enter globally into the asymptotic variance while in Theorems 2 
and [3] only the values of the weights functions at their endpoints appear in 
the asymptotic variance. 

The asymptotic variance depends on the maturity. For positive volatility 



this dependence is through d in (3.4). The martingale condition is equivalent 
to (7^/2 + 7 — A + e^i^{x)dx = 0, especially it holds that A — 7 — cr^/2 > 
with equality if and only if A = that is in the Black-Scholes case. In the 
case of positive volatility a the asymptotic variance grows exponentially as 
r — )■ 00 if the jump intensity A is positive and grows quadratic as T — t- 0. 

For wl{l), w]{l), wl{l), wl,{l) £ 1R\{0} Theorem[2]describes the asymp- 
totic distribution of the leading stochastic error term of ct^, 7, A and fi{x), 
X ^ 0, i.e., the other stochastic error terms are of smaller order. The 
variances in Theorems [2] and [s] converge by (3.5) and by the definition of 



Ze,u{x)- If one only considers the stochastic errors of a , 7, A and /i(0), 
then the covariances converge, too. But for x 7^ the covariance of the 
stochastic errors of fi{x) and of a'^ does not converge. The same holds for 
the covariance of the stochastic errors of fi{x) and 7 as well as and A. 
The phenomenon that the covariances do not convergence comes from the 
fact that the stochastic error centers more and more at the cut-off frequency. 
The sequence of cut-off values has a crucial influence on the covariance. For 
estimators of the generalized distribution function of the Levy density this 
is likely to lead to a similar dependence on the sequence of cut-off values as 
observed in [24] for deconvolution with a supersmooth error density. 

4 Applications 

4.1 Construction of confidence intervals and confidence sets 

For (T = we define confidence intervals 



[A(0) - s^(Q)eU^/'^qa/2, /i(0) + s^(o) 
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where x G ]R\{0}, qa denotes the (1 — a)-quantile of the standard normal 
distribution and 



'7 



1/2 
1/2 
1/2 



5(0)(/o^'^«o(n)'d7x) /(27r) 



with s{x) := l^pTx^ix^T^f) exp(T(A — 7))/T. We fix some arbitrarily slowly 
decreasing function h with h(u) — )• as |n| — oo. We denote by HsiR, fmax) 
the subset of Levy triplets in GsiR, o"max) that satisfy in addition 



II-^mIIoo < -R, \Tn{u)\ < Rh{u), VmER. 



(4.2) 



The additional conditions are used to extend the convergence in the theorems 
to uniform convergence, see Theorem |4j 

Corollary. Let a = 0. On the assumptions of Theorem^ and on the as- 
sumption that 6 is positive and continuous 



lim inf 



a 



holds for the intervals (4.1) and for all p G {j, X, fi{x)\x G R}. 

Remark. We can take the two parameters 7 and A and define := {(7 + 
eU^^'^s^x, X-\-eU^^'^sxy)'^\x'^+y'^ < ka}, where ka denotes the (1— a)-quantile 
of the chi-square distribution X2 with two degrees of freedom. Then it holds 



lim inf 



e ^e) = 1 



a. 



For cj > we define confidence intervals 



^M(0),e 



[7 - S.eU-'e^'^'^y^q^/,, 
s^eue I qa/2, 
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I s .rr-l„TCT2c/2/2 ] 

+ s^eu e (la 2]^ 



[A - sxee^- ^ '^qa/2^ 
[(i{x) 



X + sxeeT-'uy2g^ 



■J2\, 

■2tt2 I 



A(0) + s^(0)eC/e'^-'^V2^^/2], 
Kx)+s,eUe'^'^'u'l^q^/2l 



(4.3) 



where x G 1R\{0}, 



L2(]R) 



exp(r(cjV2 + 7- A))rV2 
and qa denotes the (1 — a)-quantile of the standard normal distribution. 



ko(l)|/(2vr) 
V ki(l)|/(27r) / 



11 



Corollary. Let o" > 0. On the assumptions of Theorem^ 



lim inf Qr (p G /„ J = 1 - a 

e-S>0 T 



/loWs for the intervals (4.3) and for all p E {7, A,/u(x)|x G R}, where the 
infimum is over all T G 'Ws(i?, o"max) volatility a. 

For (7, A)^ a uniform confidence set may be obtained similarly as in the 
case a = 0. Since for x G 1R\{0} the covariance of and V^^ij(^^) 

and the covariance of ^^ [/(^^(x) and VF^^(7(e) do not converge, confidence 
sets for {'y,fi{x))~^ and {X, p{x))~^ have to be constructed differently. Let us 
illustrate how to proceed in this case by constructing a confidence set for 
{lJ.{xi), p{x2))^ , xi,X2 G ]R\{0}. By Theorem [s] the convergence 





holds for e 

cos([/xi) sm{Uxi) 
cos(C/x2) sin(C/j;2) 

and observe that the components of Mj^^ are bounded for U for which the 
absolute value of the determinant is bounded from below by some c > 0, 
i.e., I sin(L''(a;2 — xi))| > c. For such U{e) 

'^uU ( AfLjx^) \ dwlil) ( W,^u^,^{xi) 



^^Ta'^U{.Yl2 y [/(g) AA(X2) j 2tT \ V,^u^,){x2) 

holds for e — )• 0. By Slutsky's lemma we obtain 

1 f A/i(xi) >j ( W 



s^eU{e)eT-'U{e)y2^'^U{e) ^ A/i(x2) J \V 
for e — )• such that | sin(C/(e)(x2 — > c. We define 

r — ( ^(^1) \ ^( cos([/(e)xi) sin([/(e)xi) \ 

■ \ fi{x2) C0S(C/(6)X2) Sin([/(6)X2) ^ 

where Bf, := {s^eC/(e)e^°'^^('^)^/^(x, y)^| x^ + < ka} and /cq, denotes the 
(1— Q!)-quantile of the chi-square distribution x| with two degrees of freedom. 
Then 

, ■ .... ^riiKxi), l^ix2)V G C7,) = 1 - a 

\sm{U{€){x2-xi))\>c 

holds for ah T G GsiR, CTmax) n {fi > 0}. 
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4.2 Inference on the volatility 

The results on asymptotic normality allow to test for the value of the volatil- 
ity. Let do G [OjCTmax] and define the hypotheses Hq : T G ^s(-R, o"max) H 
{(J = (To} and the alternative Hi : T G T-LsiR, <^ma.x) H {\a — ao\ > r} for 
some r > 0. For a G (0, 1/2] there is a test, which reaches asymptotically 
the level a, i.e., lim£_^.o sup-7-IE7-[(^o-o] = ck> where the supremum is over all 
T in the hypotheses Hq. Moreover, the error of the second kind vanishes 
asymptotically, i.e., lime_s.o sup7-E7-[l — 99^0] = 0, where the supremum is 
over all T in the alternative Hi. This family of tests can be used to con- 
struct a confidence set for a. The precise construction of the test and of the 
confidence set is given in Section [7j 



4.3 A numerical example 

We consider the Merton jump diffusion model [18], where the jump density 
is specified by 

A / {x-r]y\ 

with parameters a,X>0,v>0,r]GWi and where 7 G H is determined by 
the martingale condition. We simulate data with the parameters a = 0.1, 
X = 5, r] = —0.1, V = 0.2, which implies 7 = 0.379. The interest rate 
is taken to be r = 0.06. We observe prices of n = 100 European options 
with maturity T = 0.25. The strike prices are obtained from sampling the 
data points (xj) from a centered normal distribution with variance 1/2, so 
that more strike prices are sampled at the money than in or out of the 
money. The observation error is chosen to be a centered normal distribution 
with variance 6'^0{xj)'^, 6 = 0.01. Belomestny and ReiB ^ describe the 
implementation of the estimation method in detail. 

We interpolate the corresponding European call prices linearly. The 
weight functions are chosen as in [3] with smoothness parameter s = 2, 
there denoted by r. In the simulations the confidence intervals based on 
the asymptotic distribution turn out to be to conservative. The asymptotic 
confidence intervals are based on the asymptotic variance of the linearized 



stochastic errors that is the stochastic errors, where the linearization (3.2) is 



used. Instead of taking the asymptotic variance of the linearized stochastic 



errors we derive confidence intervals from the finite sample variance (6.17) of 
the linearized stochastic errors. In the finite sample variance we substitute a, 
7, A, ^ and u by their respective estimators. This yields feasible confidence 
intervals. We bound the oracle cut-off value at 30 and perform 1000 Monte 
Carlo iterations. The coverage probabilities of 95% confidence intervals for 
<T^, 7 and A are approximately 0.99, 0.95 and 0.97, respectively. We see that 
the that the coverage probabilities are close to the prescribed confidence 
level. 
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-1.0 -0.5 0.0 0.5 1.0 

X 

Figure 1: True Levy density (blue, solid) with pointwise confidence intervals 
(red, dashed) and 100 estimated Levy densities (grey). 



Figure [T] illustrates the true Levy density and the pointwise 95% confi- 
dence intervals based on the finite sample variance. At zero there is a nega- 
tive bias since the peak is smoothed out. Furthermore, the theory predicts a 
larger confidence interval at zero involving the effective weight function wq. 
This phenomenon extends in the simulations to a neighborhood of zero. The 
precise construction of the confidence intervals and more simulation results 
are contained in [20j. 

5 Conclusion 

We have shown asymptotic normality in a nonparametric calibration method 
for exponential Levy models. These results were used to derive confidence 
intervals and confidence sets as well as to construct a test on the value of the 
volatility. We have seen in a numerical example that confidence intervals 
based on finite sample variance perform well in terms of coverage probabili- 
ties. The confidence intervals extend the calibration method beyond a pure 
point estimate and enable an assessment of the calibration error. Although 
parametric models might be fitted better, the parametric approach is always 
exposed to the risk of model misspecification and the obtained confidence 
results should be used for a goodness-of-fit test. 

The estimation method and the asymptotic normality results may be 
adapted to other models as long as there is an equation relating the option 
function to the characteristic function and the parameters of interest appear 
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in the characteristic function. The constructed confidence intervals and sets 
may be used to quantify the errors in pricing, hedging and risk management. 



6 Proof of the asymptotic normality 

First, we write A7, AA and A/x(j;) similarly as in ( |2.8[ ) for A(5"^: 

A7 = -Act^ + 2[/-i /" Im{J^fi{Uu))w}^{u)du (6.1) 
Jo 

+ 2U-^ [ lm{Aij,{Uu))w}^{u)du, 
Jo 

AA = A(tV2 + A7-2 / Re{J='i^{Uu))wl{u)du (6.2) 

Jo 



Afi{x) = UT-^ [AiIj,{Uu)wI{u)\ (Ux) (6.3) 

+ ^UT-' [{Uu-ifwliu)] (Ux) 

- iA^UT-^ [{Uu - i)wl{u)] (Ux) 
+ AXUF-^ [wl{u)\ {Ux) 

- UF-^ [(1 - wl{u))Fii{Uu)] {Ux). 



In (6.1) we can substitute Act using (2.8) and obtain two error terms in- 



volving F^ and two error terms involving A^e- By similar substitutions in 



(6.2) and (6.3) we see that all error terms either involve F^ or A^/j^, which 
we will call approximation errors and stochastic errors, respectively. 

We will now state the conditions more precisely on which the regression 



model (2.1) and the Gaussian white noise model (3.1) are equivalent. We 
restrict the Gaussian white noise model to a growing sequence of intervals 
[an, /3n] and assume as a simplification that the observations in the regression 
model are equidistant on these intervals with mesh size A„. We assume 
that {j3n — an)An — )• 0. We suppose 5^ > to be an absolutely continuous 
function and -g^^og5{x) < C to hold for some C < 00. The functions O 
are uniformly bounded 0{x) = S~^C{x,T) — (1 — e^)"*" < 1 and uniformly 
Lipschitz \0'{x)\ = 1/!^ 0"(x)dx-l|,>o} + e(^-^^'^l{x>7T,.=o}l < 4 + 
where we used Proposition 2.1 in [2] and I7I < R. These properties of O 
are used to apply Corollary 4.2 in [3]. Then e corresponds to (/?„ — a^)! ^/n, 
especially for logarithmically growing intervals one loses only a logarithmic 
factor. 
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6.1 Proof of Proposition [T] 

First we define X{u) := f^^e'^^^6{x)dW{x). We assumed that there is an 
p > 1 such that + \x\)^6{x)'^dx < oo. It is shown in |19j that there 

exists a number c > such that Y^E[|X(ti) — < c|m — u |° for all 

u,v € Wi with a := min(p/2, 1) G (1/2, 1]. Denote by Np{I,r) the covering 
number, that is the minimum number of closed balls of radius r in the 
metric p with centers in / that cover /. We define p{u,v) := c\u — vl" and 
d{u, v) := Y^E[|X(u) — A ball of radius r in the metric p covers an 

interval of length 2{r/c)^/". Thus, it holds 

Np{[-U,U],r)=\u{c/rf''' 

The radius of the smallest ball with center in [—U, U] that contains [—U,U] is 
cU°' with respect to the metric p. There exists D < oo such that d{u, v) < D 
for all M, u G R. For U large enough such that cf7" > D we have the entropy 
bound 

/•oo 

J{[-U,U],d)= ilog{Na{[-U,U],r))f^dr (6.4) 
Jo 



D 



{log{Nd{[-U,U],r))f^dr 



D 



< / {logiN,i[-U,U],r))f'dr 



< (log(2C/(c/r)^/"))^^^dr 

< (log((2"C/"c/r)i/")y^^dr 

< oT^I^ / (log(2"?7"c/r))^/2^r. 



here we substitute r = 2^U'^cs, 

pD/{2°'U"c) 

< a-^/22«[/"c / (log(l/s))^/2ds. (6.5) 

Jo 

This integral is solved by 



Vlogy-idy = ^ — Erf(\/logx-i) + x\J\ogx~^, 

where Erf(y) = e^^^ dt. For all y > the estimate 1 — Eif{y) < 

exp(— y^)/(y^y) holds. For each a > this yields c > such that for all 
X G]0,l/2°[ 

/•X 

/ \/log y^^dy < cxy/logx~^ . 
Jo 
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Thus, (6.5) can be bounded by 



{a-^/^cD).J\og{2^U^c/D) < ^log{U) 



as U — )• oo. Thus, y^log{U) is an asymptotic upper bound of the entropy 
integral (6.4). We apply Dudley's theorem pLS, p. 219] to the real part of 
X. For all g > 1 this yields a continuous version X' of Re(X) with 



E 



sup |X'(n)|'? 

uG[-U,U] 



< (log(C/))^/' 



(6.6) 



as [/ — )• oo. Since X' and Re(X) are both continuous they are indistin- 



guishable and (6.6) holds for Re(X) likewise. We obtain analogously for all 



E 



as U — >• oo. We estimate from above for all g > 1 



sup |Im(X(u))|« 

ue[-u,u] 



E 



sup \Ce,u{uW 
u£[-l,l] 



< sup 

uGl-U,U] 

< 



eiu{l + iu) 



T{l + iu{l + iu)TO{u)) 
eUVl + C/2 



E 



sup \X{u)\'i 

uel-U,U] 

\ 

E 



rexp(r(-a2C/2/2 + + 7 - A - ||-F^||oo)) 



sup \X{u)\'' 



< (eC/Vlog(?7)exp(ra2c/V2) 
as [/ — )• 00. 



(6.7) 



6.2 The linearized stochastic errors 

The linearized stochastic errors are of the form fj{u)Ce^u{u)du, where 
fj with j = l,...,n are Riemann-integrable function in L°°([0,1]). Next 
we will show that these are jointly normal distributed. Almost surely C^^jj 
is continuous. Thus, almost surely the ^jCe^u are Riemann-integrable and 
almost surely 

as m — >■ 00. Let C > be such that \\fj ||oo < C for all j = 1, . . . , n. For each 
m the n sums are joint, centered normal random variables. For m — )• 00 the 
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covariance matrix converges by the dominated convergence theorem with 
the dominating function C^sup„g[Q ;^] \C^^jj{u)\'^, where sup„g[o,i] \^fi,u{u)\'^ 
is an integrable random variable by Proposition [T] Thus, the characteristic 
function converges pointwise. By Levy's continuity theorem this shows that 
the sums convergence jointly in distribution to normal random variables. So 
fj{'^)^e,u{u)du are jointly normal distributed. 

For a fixed cut-off value U the linearized stochastic errors are jointly 
normal distributed. So the natural question is whether the appropriately 
scaled covariance matrix converges for [/ — )• oo. 

Let Wj,Wk ■ [0, 1] — 7- IR be Riemann- integrable functions in L°°{[0, 1]). It 
holds 



-1 

T I u;j(u)£,,c/(n)e-^^"^^dn 







Jo J-oa 

where 

f . Wj{u){-u'^ + iu/U) 
fu[U) ■■= frr-r- (TT w (^-^^ 

exp(ry^/i((/ti)j 

and 9j := Xj + Ta'^ + T7. We define analogously 

Wk{u){-u'^ + iu/U) 

9u{u) := (6.9) 

exp(i J- ii{Uu)) 

and Ok := Xk + Ta"^ + T'^. We extend fu and gjj by zero outside the interval 
[0,1]. 



Since E 



< 00 we may apply Fubini's theorem and 



then we apply the Ito isometry to obtain 

j.2g2T(<72/2+7-A)jg 



1 rl 

Wj{u)Ce,u{uy^''''"^u / ^^;fc(^;)£e,c/(^')e-*^''^■*d^; 



J 
e^U^ / / / /^(,x)e*^"(^-^^)+^'^'^'"'/2 

JO Jo J -00 



^^(^)eit/»;(x-efc)+Ta2c/2^2/25(3;)2dxd'udv. (6.10) 



To separate real and imaginary part we will also need 
-1 /•! 



2.2g2T((72/2+7-A)jg 



( u;j(u)£,,c/(u)e-^^"^^dn f Wk{v)C,^u{v)e-'^''''^dv 
' Jo ' 

e^U^ / / / /t;(n)e^^"("-^^)+^^'^'"'/2 

Jo Jo J-00 

gu{v)e'''''^^-'^^+^'^'''''"'/^6{xfdxdudv. (6.11) 
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Lemma 1. Let o" = 0. For j = 1, . . . ,n let Xj £ IR and let Wj : [0, 1] — )• R 6e 

Riemann-integrable functions in L°°([0, 1]). Let 6 be continuous at xi + Tj, 
X2 + Tj,...,Xn + T-f and let e L^{WL). Let . . . ,W^^,V^„ . . . ,V^„ 
be Brownian motions. If Xj = Xk let Wx^ = Wxj, and Vx^ = Vx,. otherwise let 
the Brownian motions be distinct. Let the set {Wx^ , • • • , Wx„ , Vx^ , ■ ■ ■ , 
consist of independent Brownian motions. Then 

converge jointly in distribution to 



rexp(r(7-A)) vJo JO 
as U — )■ oo. 

Proof. We will first consider the We have seen that 



2.2g2T(7-A)jg 



1 pi 

Jo 



e^U^ I Tfu{U{x - ej))Tgu{U{x - e,))6{xfdx, 



where fjj and gu are defined as in (6.8) and (6.9), respectively, and 



Tfu{y)Fgu{y)KylU + ejYdy. 

-oo 



We notice that Fd"^ G L-'^(]R) impHes 5"^ € L°°(]R) and we obtain by the 
Plancherel identity 



2^e^U^ fu{u){gu{v)*F-\5{ylU + ejY){v)){u)du, (6.12) 



since the support of fu is [0, 1]. Because we are only interested in the limit 
?7 — >• oo we may assume U > 1. By the Riemann-Lebesgue lemma J-fi{u) 
tends to zero as \u\ — oo. The factor fu{u) converges for each u £ [0,1] 
to —u'^Wj{u) as f/ — )• oo and the functions are dominated by a constant 
independent of U. In order to apply dominated convergence it suffices that 
the second factor is dominated by a constant independent of U and converges 
stochastically with respect to the Lebesgue measure on R. 

{9uiv)*:F-\Siy/U + 9jf){v)){u) 

gu{u - v)T-H6{y + ejf){Uv)Udv 
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By assumption T^'^ lies in L^(]R) and so does T ^{5{y+6j)'^). A dominating 
constant is \^\\wk\\oo exp{T\\I'fi\\oo) W^^^i^iu + Oj)'^)\\L^(]R.)- It holds 



T-\5iy + 9jf){Uv)Udv 



T-\6{y + ejf){v)dv 
TT-\6{y + 9,fm = 6{9jf. (6.13) 



6u{v) := F^^{5{y + 9jY){Uv)U is the multiple of what is called approximate 
identity or nascent delta function. The basic theorem on approximate iden- 
tities states that h * 6n converges to h in L^(IR) as n — )• oo for h G L^(IR). 
Thus, {—v'^Wk{v)l[Q^ij{v))*6u{v){u) converges to —u^Wk{u)l^Q ij{u)6{9j)'^ for 
[/ — )• oo in L^(IR) [131 P- 28] and in particular stochastically. If w 7^ 0, then 
there is a neighborhood of u where gu{u) + ^^^'W^A:(^)l[o,i](^^) converges uni- 
formly to zero. {gu{v) +v'^Wk{v)l[Q^i]{v)) *6u{v){u) converges almost surely 
and in particular stochastically to zero. Therefore, gu{i')*^u{v){u) converges 
to —u'^Wk{u)l[(j^i]{u)6{9j)'^ stochastically with respect to the Lebesgue mea- 
sure on R. 

We obtain under the limit [/ — ?• 00 by the dominated convergence theo- 
rem 



1 



lim 



E 



2tt 5(9 j) 



T2 exp(2T(7- A)) 
T2 exp(2T(7- A)) 



-u^Wj{u)\Q^i^{u)) (-M^u;fc(u)l[o,i](u)) du 



■00 
1 



u Wj{u)wk{u)du. 



(6.14) 



Without taking the complex conjugate in (|6.12|) we obtain 
= 27r 



Wj{u)Ce,u{u)e ^^"''Jdn / Wk{v)C^^u{v)e 



-iUvx 



dv 



fu{u)igu{-v) * F-\5{y/U + 9,Y){v)){u)du. 



The same argumentation as before leads to 



Mm E 



Wj{u)C,^u{u)e-'^''='^du / Wk{v)C,^u{v)e-'^'"'^dv 



2-n5{9j 



r2 exp(2r(7- A)) 
0. 



(-u^ii;j(u)l[o,i](n)) (-u^Wfc(-n)l[o,i](-u)) du 

(6.15) 
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We combine (6.14) and (6.15) to obtain 



lim 


Re / 
Jo 


lim 

U^oo €''11 


Im / 

Jo 


7r5(%)2 




r2 exp(2r(7 


-A)) 



exp{iU uxj) 

^ Wj{u)Ce,u{u] 

Iq exp{iUuxj) 



du Re 
du Im 



^ Wk{u)Ce,u{u) 
exp{iUuxj) 

^ Wk{u)Ce,u{u) 

exp(iC/nxj) 



dv 
dv 



u'^Wj{u)wk{u)du. 



From (6.14) and (6.15) it also follows that the covariance between real and 
imaginary part vanishes asymptotically. 

In the Xk we have to show that the covariance vanishes 

asymptotically. Without loss of generality we assume xj < Xk- We define 
e := {Oj + ek)/2. Then Oj < 9 < Ok. 



2^2g2T(7-A) 



-E 



u;j(u)£e,{/(n)e-^^"^^du / Wkiv)C,^u{v)e~iU^^kdv 



Ffu{U{x - ej))Fgu{U{x - ek))5{xYUdx 



/oo 
^fu{y + U{e - 9j))Tguiy + Ui9-ekmy/U + 9fdy 
-oo 

By the Plancherel identity and by the dominated convergence theorem 

J^fu — T {^—u^Wj{uj) and Fgij — F [—u^Wk{u)) 

in L2(1R) for [/ — )• oo and especially the ^^(IR) norms converge. From the 
assumption FS"^ G L^(IR) follows that (5^ G L°°(1R). By the Cauchy-Schwarz 
inequality 



lim 



/ J'fuiy + U{e - 9j))Fgu{y + U{9 - 9k))5{y/U + 9fdy 



< lim 

U^oo 



I oo 



U\\l'2(!R.) 



uie-ek) 



1/2 



\J^guiy)\ dy 



0. 



A similar calculation shows that the integral over (0, oo) converges to zero 
and consequently 



2.2g2T(7-A) 

lim 277^ ^ 

u-^oo e^U-^ 

The same way follows 

2.2g2T(7-A) 

hm ^— ^ E 



r "dZi 

exp{iUuxj) Jq exp{iUvxk) 



-dv 



0. 



^ Wj{u)Ce,u{u) Wk{v)Ce,u{v) 

Jq ex.p{iUuxj) Jq ex.p{iUvxk) 



0. 
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The l/{eU^/'^) 'Wj{u)C^^u{u)e~'^^'^^idu are centered normal random 
variables and their covariance matrix converges to the covariance matrix 
of the claimed limit. Thus, the characteristic function converges pointwise. 
By Levy's continuity theorem this shows the convergence in distribution. □ 

Lemma 2. Let a > Q and 6 G L~(1R). Let wu,wu G L°°([0, 1],C) be 
Riemann-integrable and let there be a constant C > such that for allU > 1 
\\wu\\oo, WwuWoo ^ C- Let there be a,d : [l,oo) — t- C such that the condition 

lim sup sup \wu{u) — a{U)\ = (6.16) 
U>lu£[l-5/U,l] 

and the corresponding condition for wu and d hold. Then 



J^(e2 exp(r(j2C72) 
a 



E 



E 



0, 



{UyaiU) JZ, diyfdy 



0. 



exp(2T(a2/2 + 7 - X))T^a'^ ^ 

Remark. Obviously a{U) := is the only possible definition. Thus, a 

describes the dependence of wu on U at one. 



Proof. We notice that (6.10) applies to the complex- valued functions and 
yields for Wj := wu and := wu with the definitions (6.8) and (6.9) of fjj 
and gjj, respectively, and with 6 := Ta"^ + T-y 



2.2g2T(<72/2+7-A)]g 



^{fuiu)e 



wuiu)Ce,u{u)du / ■wu{v)Ce,u{v)dv 



){u{x-e)) 



fOO 



00 

■2rr2„2 , 



gu{v)e^'' ^ ^ /2 * F-\5{y/U + eY){v)^ {u)du 

= 27Te^U^ f [' fu{u)M^) 
Jo Jo 

T-\5{y + ef){U{u - t;))e^"'^'("'+"')/2^«d7;. 
For all (5 > we have 



(6.17) 



lim ra2[72e-^-'^V2 
I/— ^00 



ue 



hm e-^-'^'/2 

[7— >co 



1-5/(7 
1 



l-S/U 
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1 - lim e-^'^'^'(i-(i-'5/^)')/2 
1 - lim e-^'^'^'(25/c/-52/c/')/2 

1 - lim e-^'^'^^+^'^''^'/2 = 1. (6.18) 

U^oo 



For the product of two such sequences we obtain for all 5 > 



lim rV^C/^e-^-'^' uve^^'^'^^'+^'y^dudv = 1. (6.19) 

Jl-5/U Jl-5/U 



Likewise 

lim T^a^U^e-^'^'^' C uve^^'^^'^^'+^'^'^dudv = 1 (6.20) 



/O JO 

holds. We scale the integral in (6.17) appropriately: 



lim (r^a^U'e-^'^'^' [ [ fu{u)^) 

'-^<x^\ Jo Jo 

T-\5{y + eY){U{u - i;))e^"'^'("'+''')/2^nd?; 

-a{U)~a{U)—j 5{yfdy\ (6.21) 



^0 



F-\5{y + 0)2)(C/(7X - ?;))/c/(^i)5t/(t^)/(™)d7xdi; 



1 



-a{Uya{U)—j 5{yYdyj (6.22) 

We recall that in the Gaussian white noise model we assumed 5 to be 
in L2(II). Since F-^{5{y + eY){U{u - v))fu{u)gu{v)/{uv) is bounded in 
L°°([0, 1]2) independently of C/ for [/ > 1 and since the difference between 
( |6.20 ) and ( |6.19 ) is zero, only the integral over [1 — 5/U, 1]^ contributes to 
the limit. For all 5 > the limit equals 



Jl-S/U Jl~5/U 



:F-^5{y + e)^){U{u - v))fu{u)gu{v)/{uv)dvdu 
-a{U)~a{U)—j 5{yfdy\=Q, (6.23) 

which can be seen the following way. F~^[5{y + 9)"^) is continuous and 
\U{u — v)\ < 6 for all u,v £ [1 — 6/U, 1]. So by choosing 5 small enough 
J^-^i5{y + 9)'^){U{u - v)) gets arbitrarily close to T-^i6{y + 61)2) (0) = 
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{1/2tt) J°°^6{y)'^dy. By (6.16), wu{u) tends to a{U) and wu{v) tends 
to a{U) for 5 tending to zero. By choosing U large the factor {—u + 
i/U)/ exp(TJ^/i(?7n)) gets close to minus one for all m G [1 — (5/C/, 1]. Thus, 
for small 5 and large U the term fjj{u)gu{v)/{uv) is close to a{U)a{U) for 
ah G [1 - 5/CAl]. 

Rescaling ( |6.17 ) and taking the limit f7 — t- oo leads to 



lim . 

c/^ooVe2exp(ra2C/2) 



E 



a{Ura{U) r^5{y?dy 
exp(2r(cj2 /2 + 7 - A))r4(j4 
/ 27r;74exp(-rc72C/2) 



1 /.I 




fu{u)gu{v) 



lim ( — — 

l/^ooVr2exp(2r(o-2/2 + 7- A)) ^0 

27r 



a(?7)a(C/) 
exp(2r(c72 /2 + 7 - A)) rV427r 7_ 



5(y)^dy = 0, 



(6.24) 



where we used that (6.21) is zero. By (6.11) we have 



Hfu{u)e^''^'-'"){U{x-e)) 



Hfu{u)e'''''-'''){U{x-e)) 



T{gu{-v)e'^"^^^^^/'^){U{x - 6))5{xfdx 



fu{u)e 



(gu{-v)e^"^^^''^''^ * F-^{5{y/U + ef){v)^ {u)du 

27re^U^ j^' j\u{u)gu{-v) 

F-\5{y + ef){U{u - ?;))e^"'^'("'+''')/2^7;dn 
1 rl 
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F-\6{y + ef){-U{u + v))fu{u)gu{v)/{uv)dvdu. 




uve 



T(t2C/2(m2+i,2)/2 
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Rescaling as in (6.24) leads to 
lim ^— ttK 



2ttU^ exp(-rcj2[/2) 



lim 



1 fi 



^0 



uve 



u^oo r2 exp(2T(f72 /2 + 7 - A)) 
T-\6{y + ef){-U{u + t;))/t7(n)5c/(i;)/Mdnd7; 

since T^^{6{y + 9f'){u) — ^ for \u\ 







Tf72(72(«2+t.2)/2 



oo. 



Lemma 3. Let o" > anrf 5 G L°°(1R). Let xq £ H and for j = 1, . 

Wj : [0, 1] — >• R 6e continuous at one, Riemann-integrable and in 
Then 

1 

' u;j(n)£,,[/We-^^"^''dn 



(6.25) 

□ 

. , n Zet 
10,1]). 



eexp(r(T2[/2/2) 
converge jointly in distribution to 

¥\\L'^{^)Wj{l) 



^/2exp(r(c72/2 + 7- A))r2cj2 

asU — )• oo, where W and V are independent standard normal random vari- 
ables. 

Proof. We define n;(7(n) := Wj{u)/ exp(iC/nxo), wu{u) := Wk{u)/ exp(iC/nxo), 
a{U) := Wj{l)/ exp{iUxQ) and a(C/) := Wk{l)/ exp(iC/xo) and apply Lemma [2] 
Condition (6.16) is satisfied since Wj and vuk are continuous at one and 
since exp(— if/nxo) = exp(— iC/xo) exp(if7(l — u)xo) where U{1 — u) < 5 
for n G [1 - 5/U,l]. We note that a{U)a{U) = Wj{l)wk{l) is real. By 
Lemma [2] the covariances converge to the covariances of the claimed limit. 
The convergence in distribution follows by Levy's continuity theorem. □ 

Lemma 4. Let o" > and 5 G L°°(1R). Let wi,W2 ■ [0, 1] — )• IR 6e Riemann- 
integrable, in L°°([0, 1]) and continuous at one. Let xi,X2 G H and denote 
X2 — xi by if. Then 



wi{l) 



1 

as U — )■ 00. 
Proof. We define 



^ W2{u)Ce,u{u) 



dn 



•^2(1) 



^ wi{u)C^^u{u) 



r,iU uxi 



dn 



0, 



VJl{l)w2{u) 



W2{l)wi{u) 



exp{iUux2) exp{iUip) exp{iUuxi) 
wjj{u) fulfills condition ( 6.16| ) with a{U) = for all f7 > 1. Lemma [2] yields 



lim 



-E 



wif{u)Ce,u{u)du 



and the statement follows by Levy's continuity theorem. 



□ 
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6.3 The remainder term 



In this section, we show that the contribution of the remainder term to the 
estimation vanishes asymptotically. We recall that the remainder term TZ^^u 
depends on the Levy triplet. 

Lemma 5. Let do > 0. Let wjj S L°°([0,1],C) be Riemann-integrable and 
let there be a constant C > such that \\wu\\oo < C for all U > 1. If 
eC/(e) Vlog(f/(e)) exp{Ta'^U{e)'^ /2) -^0 ase^O, then for all Levy triplets 
with o" < o"o 



1 



'^;7{e)(^)'^e,C/(e)(^)d^i 0, OS e — )■ 0. 



eexp(r^72c/(e)2/2) ^ 

Proof By the identity 7^e,c/(u) = (l/T) log>^[7(„)(l + TC^^u{u)) - Ce,u{u), 
where kF [u) < 1/2, we have to show that for U = U{e) 

\ / wu{u) (log>«^[/(„)(l + T£.,^u{u)) - TC,^u{u)) du 

(6.26) 

converges in probability to zero. 

For 2 G C holds log(l + z) — z = 0(|zp) as \z\ — > 0. We define g by 
g{z) := (log(l + z) — z)/\z\'^ for z ^ and ^(0) := 0. There are M and 
T] > such that \g{z)\ < M for all \z\ < r]. We may assume that rj < 1/2. 
If the logarithm in the definition of g is replaced by the trimmed logarithm 
log>^ with K £ (0,1/2] then g remains unchanged for \z\ < 1/2. Thus, 
the statement holds uniformly for all ^^(-z) := (log>K(l + -2) — ^)/\^\'^ with 
kG]0,1/2]. 

By Proposition [1] we have sup„g[_x,i] l^e,c/('u)| 0. Let r > be given. 
Eventually we have 

3u G [-1, 1] : I log>,a(,)(l + TC,,uin)) - TC,,u{u)\ > MT^\C,,uint 

< Q I T sup \Ce,u{u)\ >7]] <T. 

\ ue[-i,i] J 
Except on a set with probability less than r we have eventually 
1 



eexp(rag;72/2) 



'0 



j wuiu) (^log>^u(„)(l + TCe,uiu)) - TCe,u{u)^ du 



< Z\rr2/n^ I \wu {u) C,,U {uf\du. (6.27) 

eexp(rcj^C/V2) Jq 



Hence (6.26) converges in probability to zero if (6.27) converges in proba- 
bility to zero. The convergence 

1 

'wu{u)£.e,u{u)'^\du 



eexp(Ta2[/2/2) JO 
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holds even in since 
1 



< 



< 



eexp(rcjgC/2/2) 
C 

eexp(rcjgC/2/2) 
C 

eexp(rcjgC/2/2) 



ii;[7(u)£e,;7(?^)^|dn 
/ |£e,c/(w)^|du 

JO 



(6.28) 



r(l + iUu{l + iC/n)J^C'(C/n)) 



E 



iUux 



5{x)dW{x] 



du 



< 



C 



1 e\U'' + U^)uGicp{Ta'^U'^u^ 



2 

L2(]R) 



eexp(rcj2C/2/2) Jq r2 exp(2r(a2/2 + 7 - A) - 2r||^^|| 
for (7 = this converges to zero and for o" > we further calculate 
Ce(l + f/2)||5||i2(j^) /o 2ra2C/2uexp(r(T2t/2n2)d 



dn, (6.29) 



in 



exp(rcjgC/2/2)2rV2 exp(2r((j2/2 + 7 - A) - 2r||J^/z||c 



Ce(l + C/^ 



2,(^)(exp(Ta2[/2) 



exp(rCT2c/2/2)2r3cj2 exp(2r(a2/2 + 7 - A) - 2T\\T^JL\l 

2rr2 , 



< 



Cml,^^^e{l + U^){eMT4U'm) 



2T3cj2 exp(2r(fT2/2 + 7 - A) - 2r||J^^||c 



as e ^ 0. Thus, (6.26) converges in probability to zero. 



□ 



Lemma 6. Let wjj E L°°([0, 1],C) be Riemann-integrable and let there be 
a constant C > such that WwuWoo < C for all U > 1. Lf U{e) — )• 00 and 
eC/(e)^/^ —7-0 as e — >■ 0, then for all Levy triplets with cr = 

I Wuie){u)n,^uie){u)du ^ 0, 



as € 



0. 



Proof. We follow the proof of Lemma 5 sup„gr_]^ ii \^e,u{u)\ ^ holds by 
Proposition [l| W e set do = and divide by C/^/^ in (6.26) and (6.27). Then 
we use that (6.28) is bounded by ( 6.29| ), where we set do = a = and divide 
by C/3/2 again. We obtain 



.E 



\wu{u)Ce,u{u)'^\du 







Ce(C/i/2 + ^5/2)11^11 



L2(]R) 



eC/3/2 

as e — )• 0, which implies the desired convergence. 



^ r2exp(T(2(7-A)-2||^M||lo)) 



□ 
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6.4 The approximation errors 

The approximation error can be controlled as in [2] using the order conditions 
(2.9) on the weight functions. The Levy triplet T = ((7^,7, /i) was assumed 
to be contained in Gs{R-,(^ma.x)-, especially /i is s-times weakly differentiable 
and maxo<fc<s Ha*^^^ IIl2(r) < R, 11/^^'''' lU < R- 

We use {iuyj^fj,{u) = J^n^'^^u) and the Plancherel identity to bound the 
approximation error by 



Jo 



-1 
00 



fi'^'\x/U)U-^T-^{wl{u)/{iUuy){x)dx 
Analogously we obtain 



1 

IP 
27r 
IP 



(6.30) 



lm{Tfi{U u))w}^{u)du 



2 / Re{J='n{Uu))wl{u)du 



<[/-(^+2)||/.W||oo||-F(^^(n)K)|Ui(I,), 

(6.31) 

<U-^^+'^\\f,^^^\U:F{wl{u)/un\\LHU)- 

(6.32) 



The last error term in (6.3) can be bounded by 

\UT-^ [{l-wj,{u))Tfi{Uu)] {Ux)\ 
U 



2tt 



{I - wj,{u))Tfi{Uu)e-'^'''^du 



1 



27rC/^-i 



1 — w},(u) , s 

{e*^"^--F^(")(?7n)dn 



< 



Ir' Woo 



T 



-1 



\-w],(u) 



1-wUu) 



(6.33) 



6.5 Proof of Theorem [T] 



We consider the decompositions of the estimation errors (6.1 )-(6.3). The un- 
dersmoothing eC/(e)(^''+^)/^ — )• 00 is equivalent to C/(e)~(^+^ = o(eC/(e)~^/^) 
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and implies in view of (6.30) that the approximation error of a is asymp- 



totically negligible. Likewise the error terms 2C/ ^ lm{J^fi{Uu))w]y{u)du, 

2 £ Re{Tn{Uu))wl{u)du and UT'^ [(1 - wl,{u))Tn{Uu)] (Ux) can be bounded 
as in (6.31), (6.32) and (6.33) and vanish asymptotically. Since cj^ converges 



with a faster rate than 7 and 7 converges with a faster rate than A, the error 
Aa"^ vanishes asymptotically in (6.1) and in (6.2) as well as A7 is asymp- 



totically negligible in (6.2). For x 7^ we can apply the Riemann-Lebesgue 



lemma to the second, the third and the fourth error term in (6.3 ) and we see 
that they are of order 0(Q(eC/(e)^/^). For x = due to the symmetry of 
the third term vanishes asymptotically but the second and the fourth term 
do not. The error terms of fi{x) we have to consider are in the case x 7^ 

UT-^ [A^j,{Uu)wl{u)] {Ux) 



and in the case x 



2tt' 




iRe (AV^,(C/M)e-*^"^) dn 



F'^ [AMUu))wl{u)\ (0) 



+ 







Re{A'iP,{Uu))wl{u)duF^^ [u^wl{u)\ (0) 



- 2 / Re{Ailj,{Uu))wl{u)duF-^ [wl{u)] (0) 
1 /-i 

—2 / Re{A'ilj^{Uu))wQ{u)du. 
27r Jo 



(2.9) implies that w]j, w}^, 'w\ and w}^ are continuous and bounded, especially 



they are Riemann-integrable and in L°°([— 1, 1]). Applying Lemma [T] to the 
linearized stochastic errors and Lemma |6] to the remainder terms yields the 
theorem. 



6.6 Proof of Theorem [2] 

To see the first line we set xi = X2 = 0, wi = 1 and W2 = Wa in Lemma |4] 
and wu = in Lemma [5] The second and third line follow analogously. 
In order to derive the last line we observe 

[AA{Uu)wl,{u)] (Ux) 

2 / Re(AV',(C/u)e-^^"^)w;i(u)dn/(27r) 



and apply Lemma |4] with xi = 0, X2 
remainder term vanishes by setting wjj{u) 



X, wi = 1 and W2 = w^^ . 
w^^{u)e~^^^^ in Lemma 5 



The 
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6.7 Proof of Theorem [3] 



By the undersmoothing eU^'^^ eicpiTa'^U'^ /2) — )• oo we have = 
o(eC/~^ exp(T(7^C/^/2)) so that the approximation error of a vanishes. A 
similar reasoning apphes to the approximation errors of the other estima- 
tors. Since a converges with a faster rate than 7 and 7 with a faster rate than 
A the leading stochastic error terms are given in Theorem [2] and the conver- 
gence of the first three lines follows by this theorem. For x 7^ all stochastic 



errors in (6.3) are negligible except the first one. We obtain the convergence 
in the last line by Theorem pi We observe that F~^[uw^^{u)\{^) = 0, since 
is symmetric. For x = Q the relevant stochastic error terms are 

J-i [/\MUu)wl{u)] (0) 

+ / Re{^iJe{Uu))wl{u)duF-^ [u^wUu)] (0) 
Jo 

-2 [ Re{Aij,{Uu))wl{u)duT-^ [wUu)]iO) 
Jo 

= -'-2 / Re(AV'.(C/n))wo(n)dn. 







We apply Lemma|4]with = 2:2 = 0, = 1 and 11)2 = wq to this term. The 
remainder term converges to zero by Lemma [5j This shows the convergence 
in the next to last line. 



7 Tests and a confidence set for the volatility 

In this section, we test the hypotheses Hq : a = ao against the alternative 
Hi : |(T — cjol >T,r>0, and construct a confidence set for a. We assume 
that T G (Tmax) especially we assume that a € [0,(Tmax]- 

For o"o > the most natural test statistic is the following. In or- 
der to apply the uniform version of Theorem [3] under Hq we choose a 
cut-off value U{e) such that e[/(e)2 log(C/(e))i/2 exp(Tcj2[/(e)2/2) ^ and 
eC/(e)*+^ exp(T(TQ[/(e)^/2) 00 as e 0. Let ct^ be the estimator cor- 
responding to this cut-off value C/(e). We ensure o"^ G [Oj^max] W taking 
the maximum with zero and the minimum with cr^g^^. Likewise we ensure 
7 G [—R,R] and A G [0,-R]. We choose a weight function with wl{l) ^ 
and define 

'''' d.oeexp(ra2C/(e)2/2)' ' exp(r(aoV2 + 7 - A))r2a2 ' 

Under Hq the test statistic S^q converges uniformly in distribution to a 
standard normal random variable by the uniform version of Theorem [3j We 
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decompose 



S, 



d„,€exp{Ta^Uiey/2) 



+ ^ 



(7.1) 



We will show that for a < (Tq — r the first term converges uniformly in 
probability to zero. The approximation error contributes a term that con- 



verges deterministically to zero by the bound (6.30) and by the assumption 
eC/(e)^+i exp(T(TgC/(e) 2/2) 



oo as e 



0. The remainder term of the 



stochastic error in the first summand of (7.1) converge uniformly in prob- 



ability to zero by Subsection 8.3 and the linearized stochastic error by the 
following lemma. 

Lemma 7. Let 6 G L°^(1R). Let wu G L°°{[0, 1], C) be Riemann-integrable 
and let there be a constant C > such that WwuWoo ^ C for allU > 1. Then 
for all ri,T > 



sup Qr 

T&Hs{R,ao-T) 



1 



;exp(ra2[/2/2) 



0, 



as [/—)• oo 

Proof. It is equivalent to consider for all r' > the supremum over 7is{R, \/ ctq 



such that cj^ < (Tq — r' is satisfied for each T. By (16.17^ we have 



1 



sup 



-E 



1 .1 
wu{u)Ce,uiu)<iu / wu{v)Ce,u{v)dv 
Jo 

1 rl 



27rC/4exp(2r(A-7-aV2)) . / ,r , , 

sup ^ , 2tt2\ / / fu{u)fu{v) 

T-\6iy + e)^)iUiu - ?;))e^"'^'("'+''')/2d^dt; 



47rL/^exp(6ri?)C2 „ 

- r2exp(rr'c/2) "-^ "~ 



0, 



as [/ — >• oo, where we used WfuWoo < V2Cexp{TR). 



□ 



We have seen that for cr < ctq — t the first term in (7.1) converges to 



zero uniformly in probability. The second term is d^^ times a deterministic 
sequence converging to — oo. We note that d^Q is bounded from above and 
below. Consequently it holds lime_>o inf-^^(ij^o-o-T) Q( 5*0-0 < c) = 1 for all 
c S R. We would like to make a similar statement in the case a > ao + t. 
Unfortunately for a > ao the variance of jC^ ;7(e)(l) does not converge to zero. 
So it is not possible to find a bound like in Proposition [T| which converges to 
zero. Consequently the remainder term cannot be controlled by Lemma [5] 
We modify the test statistic in the following way. We choose a > (Tmax and 
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let U{e) be a cut-off value with eU{eflog{U{e))^/^ exp{Ta^U{e)'^ /2) 
and eU{ey^^ exp(T(j^i7(e)^/2) — > oo as e — ^ 0. We further assume that 

t/(6)2exp(-ra2C/(6)V2) ^ ^ ^ 

This can for example be ensured by choosing the cut-off values U (e) and 



U{e) according to (3.6) with a and a > a, respectively. Let a be the 



estimator of corresponding to the cut-off value U{e). We define 

S -S I ^{ena^-a^) 
eexjp[Ta''U (e)"'/2) 

"° eexp(r^2c7(e)2/2) eexp(ra2i7(e)2/2) ' ^ ' ^ 

Under Hq the statistic 5(jq can be written as in Lemma [8| The second 
term in (7.3) converges uniformly in probability to zero. Thus, under Hq 



the modified statistic S^q converges uniformly in distribution to a stan- 
dard normal random variable by Lemma [8) For a < uq — t the sec- 



ond term in (7.4) converges uniformly in probability to zero and the third 
term is deterministic sequence converging to — oo. As for So-,, it holds 
lime_^o iiif-H^(_R,,(To-r) Q('S'cto < c) = 1 for all c G M. But now we are also 
able to make a similar statement for a > + t. Since we bounded the 
estimators Sa^ cannot diverge faster than e~"'^C/(e)^ exp(— rcTot/(e)^/2). For 



(7 G [0 

) CTmax] the second term in (|7.4l) converges uniformly in probability to 



zero. Owing to (7.2) the third term in (7.4) tends to infinity faster than the 
bound of 5cto- It holds lime_j.o inf7"Q('S'o-() > c) = 1 for all c G M, where the 
infimum is over all T G HsiR, o"max) with a > ao + t. 

Let Qa denote the (1 — a)~quantile of the standard normal distribution. 
For o"o G (0, (Tniax) we define the tests 

/ 0, if 1 1 < qa/2 / 0, if 5o,„,, > 

1 1, ii\S^,\>qa/2 ' 11, if5o_<gi-a ■ 

For (Jo = we would like to apply Theorem [TJ To this end we choose the 
cut-off value U{e) such that eU{efl^ and €J7(e)(2^+5)/2 ^ oo. Let a'^ be 
the estimator corresponding to this cut-off value U[e), where cj2 g [0,c7^ax] 
is ensured. We assume that 5 is positive on [—TR, TR] and define 



doe eexp(rcr2^^^[/(e)2/2)' 

^°- %xpTO-a))t U^"-^^^^^J • 

Under Hq : = the test statistic 5*0 converges uniformly in distribution 
to a standard normal random variable by a similar argument as for S^q- 
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We observe that the first term of 5o is nonnegative. Under Hi : a > t 



the second term of may be decomposed as in (7.4) into a part that 
converges uniformly in distribution to zero and into a deterministic sequence 
converging to infinity. It holds lime_5.o inff 0(5*0 > c) = 1 for all c G R, 
where the infimum is over all T S TisiR, (Juiax) with a > t. We define the 
test 



^0 



0, if 5*0 < Qa 

1, if > Qa 



Since we have a test on o" = do for each (Tq G [0, amax] we may use this 
family of tests to define a confidence set for a. We define := = 0} 

and obtain lime_!.o ini-j-QTi'^ ^ ^e) = 1 — a, where the infimum is over all 
T € T-Ls{R-, o"max) with volatility a. The set is not necessarily an interval. 
For (To S (0, fTmax) the cut-off value U{e) may be chosen as a continuous 



function of fio by (3.6). The estimators , 7 and A depend continuously 



on the cut-off value {7(e), which can be seen by substituting v = u/U in 



(2.3), (2.4) and ( |2.5| ) and applying the continuity theorem on parameter 
dependent integrals. Thus, S : (0, cimax) — H, cr i—;- S'o- is continuous and 
Mg n (0, cTinax) may be written as the preimage »S'~^([— (70-/2) 9a/2]) of the 
continuous function S. 



8 Uniform convergence 

The asymptotic normality results hold for each Levy triplet T G Qs{R-, o"max)- 
The speed of convergence might depend on T. To make statements on 
confidence sets and on hypotheses tests it is useful to control the speed of 
convergence uniformly over a class of Levy triplets. We fix some arbitrarily 
slowly decreasing function h with h{u) — t- as |u| — )• 00. We will show 
uniform convergence for the class TisiR, o"max) consisting of all Levy triplets 



in ^s(i?, o"max) Satisfying the additional conditions (4.2), which we recall to 
be 

II-^mIIoo < -R, \Tniu)\ < Rh{u), VneR. 

The first condition can easily be ensured by ||^||li(]r) < R- For h{u) = \u\~^ 
the second condition can be ensured by llAt^^-* ||li(]r) < R- For each Levy 
triplet in Qs{R, o'raax) the function F^{u) tends to zero as \u\ — by the 
Riemann-Lebesgue lemma. Especially for each T G Qs{R-,(^m&x) there are h 
and i?' > such that T G HsiR', fmax)- 

In the case a > some covariances do not converge. We show uniform 
convergence for the joint distribution only in such cases, where for a > 
the covariances do converge. As it turns out it is also important that 
the covariance matrix of the limit is nondegenerated. We cover uniform 
convergence of o"^, 7, A, fi{x) and of (7, A). 
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A sequence of random variables X„ converges to X in total variation if 
sup \¥{Xn eB)- F{X £B)\^0 

B 

as n — )• oo, where the supremum is taken over aU measurable sets B. For 
sequences X^^n^ -d ^ 9, we say that they converge to X in total variation 
uniformly over 9 if 

sup sup |P(X^,„ G B) - P(X G S)| ^ 

i?6e B 

as n — )• oo. Motivated by the Portmanteau theorem we say for i9 G 

9, and X with values in some metric space that converge to X in 

distribution uniformly over 9 if for all Borel sets B with P(X G dB) = we 
have 

sup|P(X^,„ G -P(X G 5)1 ^ asn^oo. 

■dee 

Theorem 4. Let the assumptions of Theorem [i] he fulfilled and let 5 he 
continuous and positive. Each marginal convergence in Theorem^ is a uni- 
form convergence in distrihution over 'Hs{R,0) if the standard deviation is 
positive and hoth sides are divided by it. 

Let a > and let the assumptions of Theorem he fulfilled. Each 
marginal convergences in Theorem ^ is a uniform convergence in distri- 
hution over all 7" G 'H<j(-R, cTmax) with volatility a if the standard deviation 
is positive and hoth sides are divided by it. 

Remark. In both cases uniform convergence in distribution does also hold 
for 7 and A jointly and in the standard deviation on the left side 7 and A 
may be replaced by their estimators. 

The following lemma may be seen as a generalization of Slutsky's lemma 
for uniform convergence. It is the key step to show the uniform convergence 
in distribution in Theorem 01 

Lemma 8. Let X^^n, Y-&,n, G 0, n G N, and X be random vectors such that 
X^^n converge to X in total variation uniformly over 9 and sup^^g P(|l^^„| > 
(5) —7- as n —)• 00 for all 5 > 0. Let Z^^n be random variables with 
sup^gg P(|Z^^„ — 1| > 5) — )• for all 5 > 0. Then Z^^^X^^n + Y^d,ri con- 
verge to X in distribution uniformly over 9. 

Proof. For 5 > we define 

B^ := {y G 1R'^| \x - y\ < 5 for some x G B}, 
B~^ := {y G R'^jx G B for ah x with \x - y\ < 5}. 
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As and B'^ are open and closed, respectively, they are Borel sets. It 
holds n5>o^'^ = B and U5>o^"'' = B° ■ For B with ¥{X G dB) = 0, we 
have P(X G B°) = ¥{X £B)= F{X G B). Consequently 

lim F{X £ B^) = lim F{X G B'^) = F{X G B). (8.1) 

<5— s>0 (5— !>0 

Let ry > be given. For all 5 > it holds 

supP(Z^,„X^,„ + y^,„ G 

< supP(Z^,„X^,„ G S-^) + supP(|ytf,„| > 5), 
■dee i?ee 

for large n the second term is smaller than i] owing to the assumptions on 

<supP(Z^,„X^,„GS'5)+r/. 

As a single random vector X is tight meaning that for each 77 > there 
is M such that F{\X\ > M) < r]. By taking the set {x\x > M} in the 
definition of uniform convergence in total variation we obtain for n large 
enough sup^gg P(|X^^„| > M) < 27]. By considering possibly larger n we 
can also ensure sup^g^/ P(|Z^_„ — 1| > 5/M) < rj. But if |X^,n| < ^ and 
\Zd,n — 1| < 1 then \X^^nZi)^n — X^^ri\ < ^- For n large enough it holds 

supP(Z^,„X^,„ + Y^^n eB)< supP(Z^.„X^.„ eB^) + ri 

■dee -dee 

<supP(X^,„ gS2'5) + 4?7 

■dee 

< F{X G B^^) + 57? 

< F{X £B) + 6rj, (8.2) 



for S small enough by (8.1). 

Let ?? > be given. For n large enough and 6 > small enough we 
obtain similarly 

inf P(Z^,„X^,„ + eB)> inf P(Z^,„X^,„ G B'^) - r? 
■dee ' ' ' ^ee ' ' 



> inf P(X^„ G B-^^)-4r] 

■dee ' 



> F{X G B-"^^) - 577 

> F{X £B)- 67]. (8.3) 



The statement follows by combining ( |8.2[ ) and (8.3). □ 

Lemma [S] outlines how to proceed in showing uniform convergence. X^^n 
will be the leading term of the linearized stochastic error, n will be the sum 
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of the smaller stochastic errors, the remainder term and the approximation 
error and Z^^^ will be the quotient of standard deviation and estimated 
standard deviation. The approximation error is uniformly controlled over 
^^(i?, cjmax) and thus over ^s(-R, o-max), too. 

The substitution of the standard deviation by its empirical counterpart 
works as follows. We fix some x and write d instead of d{x) in Theorem [l] 
to unify the notation with Theorem [3] and to treat both simultaneously. 
For 8 continuous and positive the standard deviation depends continuously 
on 7 and A through d. 7 and A are restricted to a compact set. By the 
uniform convergence and by an upper bound of the standard deviation we 
obtain sup-7-Q7-(|Ap| > 5) — )• for all > and for p G {7, A}, where the 
supremum is over all T G o"niax) with a fixed volatility a. Since d is 

uniformly continuous in 7 and A, we obtain sup7-Q7-(|A(i| > 5) — )• for all 
(5 > 0, which gives the assumption on djd corresponding to „ in Lemmajs] 
by a lower bound on d. 

By the following lemma uniform convergence in total variation of the lin- 
earized stochastic error follows from uniform convergence in each component 
of the covariance matrix. 



Lemma 9. Let X he a normal random vector with symmetric positive def- 
inite covariance matrix A G IR.'^^'^. Let X^^n, n ^ ¥i, -d ^ 6, be normal 
random vectors with covariance matrices A^^n £ H'^^'^. If A^^n converge to 
A in each component uniformly over 9 as n —)■ 00, then X^^^ converge to X 
in total variation uniformly over 6 as n —)• 00. 



Proof. We have to show that 

exp(-(x, A^>)/2) exp(-(x, A-^x)/2) 



sup sup 

B 



B det(v27r^^,n) 



det(V2^A) 



dx 



converges to zero as n — )• 00. The determinant is a continuous function of 
the components of a matrix. For all 5 G (0, det(\/27rA)) there is G N such 
that sup^gg I det(\/27r^^^„) — det(\/27r^)| < 5 holds for all n > and 
is well-defined for all -i? G ^ if n > A^. Expression (8.4) equals 



1 

sup- 
1 

sup- 



exp(-(x, A^l^x)/2) exp(-(x, A-^x)/2) 



det(v27r^^^„) 
exp(-(x,^-ix)/2) 



det(V2^A) 



dx 



Rd I det(\/27f^)| 

det(V2^^) exp(-(x,(A-i„ 



A-^)x)/2) 



det(V27rA 



dx. 



A^^n converges to A in each component uniformly over 9. Likewise A, 
converges to A~^ in each component uniformly over 9. We now addition- 



1 
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ally require 6 E (0,Amin(^ ^)/d), where Amm(^ ^) is the smallest eigen- 
value of A~^. By going over to a possibly larger we may assume that 



all n > 



A- 



)jk\ 



< 5 for all j,k = 1, . . . ,d, for all n > N. Then for 



\{x,iA^l^-A-')x)\ 




< 



j,k=i 



12) 



where in the last step the Cauchy-Schwarz inequality is used. We see that 



(8.5) can be bounded by 

e^p{-{x,A-^x)/2) 



|det(V27rA)| 
det(\/27r^)| exp 



li/2) 



det(\/27fA)|exp(- 



li/2) 



det(V2^^)| -5 



det(V27r^)| +S 



dx. 



i.6) 



The integrand converges pointwise to zero for 5 — t- 0. For fixed 5 the func- 
tion is at the same time a dominating function, which is integrable since 
A~^ — Sdid is positive definite by the choice 5 G {0, XmmiA~^)/d). By 



the dominated convergence theorem (8.6) converges to zero and likewise 



We will show uniform convergence in each component of the covariance 
matrix. By Lemma |9] this leads to uniform convergence in total variation of 
the linearized stochastic errors provided that the covariance matrix of the 
limit is nondegenerated. To this end we assume in the case a = that 6 
is positive and that u'^Wp{u)'^du > for the weight functions w^, p G 
{(T, 7, A, 0, /i}, involved. Joint distributions may further only involve more 
than one of the estimators o"^, A or /i(0), if the covariance matrix is positive 
definite, in the case a > we assume that ^ 

and ||(5^||l2(][^) > 0. By the uniform version of Lemma 5 the remainder term 
converges uniformly to zero. By Lemma [8] each of the rescaled stochastic 
error terms 



Re{/\il),{Uu))wl{u)du, 



eexp(ra2C/2/2)«;i(l) 



eex.p{Ta^U^/2)wi^{l) Jq ' 

^] 2 / ^e{^i,,{Uu))w\{u)du, 

eexp(rcr^[/^/2)tt;j^(l) Jq 
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eexp(rfT2[/2/2)u;i(l 
converges uniformly in distribution to 

Il2(]R) 



27r 



^/2||^2| 



exp(T(a2/2 + 7- A))rV' 

where is a standard normal random variable. By the uniform versions of 
Lemma [3] and Lemma [5] we also obtain that for the second and third of the 
above stochastic error terms holds joint uniform convergence in distribution 
to independent normal variables. 

8.1 Uniform convergence in the case a = 

The convergence in distribution of the linearized stochastic errors is shown 
in Lemma [1] by the convergence of the components in the covariance matrix. 
We restrict ourselves to the case xi = ■ ■ ■ = Xn and show uniform conver- 
gence in each component of the covariance matrix. This implies uniform 
convergence in total variation by Lemma [9j We assume that S is continuous 
at all X £ [xi — TR,xi + TR]. We note that fjj, gu and 6j depend on the 



Levy triplet T- The uniform convergence of the rescaled covariances (6.12) 
will be shown by the following easy lemma. 

Lemma 10. Let U^n, h, g'&,n, g-d G ^^[0, 1],<D), n e N, I? G 6*, and M > 
such that ||/^,n||oo, HSi^Hoo < M for all n £ ¥\, for all i9 £ 9. Let 
sup^ee ll/i?,n - /i?IIli([o,i],c) "-^d sup^gg Hs'^,™ - ||li([o,i],c) as 
n ^ 00. Then 

sup / \fi),n{x)g^,n{x) - f^{x)g^{x)\dx ^ as 00. 
■dee Jo 

Proof For all -d £ 6 it holds 

,ng'd,n — f^,ngA + l/ij.nfl'is — f^g-d\ 

< M\g^^n - gi}\ + M\f^^n - 

By assumption sup^ge ll/i?,n-/^i||Li{[o,i],c) as n 00 and sup^g^ \\gi},n- 
5'i9||l1([o,i],C) — ^ as n — )• 00 and the claimed statement follows. □ 



Let us verify the assumptions of Lemma 10 for the first factor in (6.12). 
fu and v?Wj{u) will correspond to and respectively. It holds 

\fu{u) + v?Wj{u)\ 

= \u^Wj{u){\ — exp{—TJ^fi{Uu))) + iu ex.p{—TT fi{Uu)) /U\ 

< u^\wj{u)\{exp{TR{l A h{Uu))) - 1) + iuexp{TR)/U. (8.7) 
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This bound does not depend on T and converges everywhere to zero. Further 
fu{u), u^Wj{u) < \/2||tt;j||oo exp(Ti2) is a bound that does not depend on U 
nor on T- By the dominated convergence theorem 

sup Wfuiu) + '"^w'i('")llLi([o,i],c) (8-8) 

as {/ — >• oo and the conditions on the first factor in Lemma fTOl are satisfied. 



To show the assumptions of Lemma 10 on the second factor, which is 
the complex conjugate of 

iguiv)*T-HSiy/U + 9,)^){v)){u) 

gu{u - v)F~\5{y + ejf){Uv)Udv, 



we apply the following lemma. It is a uniform version of the basic theorem 
on approximate identities. 

Lemma 11. Let f,f^^n G L^{^,G), n e ¥1, e 9. Let sup^gg - 
/IIli(]R,C) ~^ /'^'^ ^ ~^ CO. Let 5^^n £ -^^^(11, C) fulfill the following proper- 
ties: 

(i) There exists c > such that ||5^,n||Li(]R,C) — c for all n G N, for all 
(ii) J^oo ^^,niy)dy = Ci) for all n E N, for all £ 0. 

(Hi) For any neighborhood V of zero we have sup^gg fyc \5^^niy)\dy — >• as 
n — )• oo. 

Then sup^gg \\6^^n * f^,n - Ci?/||Li(]R,c) as n oo. 
Proof. By the triangle inequality 

sup ||5tf,„ * f^^n - Ci»/||Ll(]R.C) 

< sup ||(5^,„ * {f{t^n - /)||li(r,c) + sup ||(5^,„ * / - Cij/lli^im <c). (8.9) 
The first term in the triangle inequality ( |8.9| ) can be bounded by 

sup * {f^^n - /)||li(IR,C) < sup ||5i(,„||il(]R,(C) - /||l1(]R,C) 

< C sup \\fi)^n - /||l1(]R,C): 

which converges to zero by assumption. To bound the second term in the 
triangle inequality (8.9) proceed as in the proof of Theorem 1.2.19 in \13\ 
p. 26]. Without loss of generality we may assume that ||/||li(]R,c) > 0- 
Continuous functions with compact support are dense in L^(1R, C). Since 
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continuous functions g with compact support are bounded we obtain by the 
dominated convergence theorem 



\g{x -y) - g{x)\dx 



as y — 7- 0. We approximate / G L^(]R, C) by a continuous function with 
compact support and see that for all (5 > there is some neighborhood V of 
zero such that 



\f{x -y)- f{x)\dx < 

oo 

It further holds 

{S&,n * f)ix) - C^fix) 

POO 

= {S^,n * f){x) - fix) / 6,^^n{y)<iy 



2c 



for all y £ V. 



(8.10) 



(fix -y)- f{x))6^^n{y)dy 

-oo 

(fix -y)- f{x))5^^n{y)dy + / {f{x - y) - /(x))(5^,„(y)dy. 
We take norms with respect to x and obtain for the first term by (8.10) 



sup 



V 



(fix -y)- fix))5^^niy)dy 



Li(]R,C) 



< sup / ||/(x - y) - /(x)||ii(]R(C)|(5tf,„(y)|dy 
•See Jv 

f s s 

< sup / —\6^^n{y)\dy < - 



(8.11) 



and for the second term 



sup 



/ {fix -y)- fix))5^^niy)dy 



Li(]R,C) 



< sup / 2||/||ii(]Rc)l'^i9,n(y)|dy < -, 

i9e Jv<= ^ 



i.l2) 



where n is taken large enough such that 
sup / \Si)^niy)\dy < 



Li(IR-.<C) 



The lemma is a consequence of (8.11) and (8.12). 



□ 
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Let us first verify the assumptions ([I]), (|n]) and (iii) for T ^{6{y+6j)'^{Uv)U, 
wliich will correspond to 6^^n- We have 



T-H6{y + er){v) 



idv 



F~\5{yf){v). 



(8.13) 



By the assumptions of Lemma [T] it holds J- 5"^ G L^(IR). The equality 

\F-\5{y + e^f){Uv)U\ = \F-\5{yf){Uv)U\ 

shows that the absolute value does not depend on T and that conditions ^ 
and (iii) are satisfied. Condition ^ is satisfied by (6.13). We have 



s\xp\\gu{u) + u ti'fcHllLicicil.c) 







(8.14) 



as in the corresponding equation ( |8.8[ ) for /[/. By extending gjj and Wk 
by zero outside [0, 1] this holds in L^R, C), too. We apply Lemma 
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to 



gu{v) and v'^Wk{v), which correspond to and / in this lemma. v^Wk{v) 
is bounded and since condition (P is satisfied we see that 5{9j)'^v'^Wk{v) is 



uniformly bounded over all Levy triplets. By Lemma 10 the convergence of 



the covariances (6.14) holds uniformly. The convergence in the analogous 



equation without conjugation (6.15) holds uniformly, too. 



8.2 Uniform convergence in the case a > 

Let us first fix o" > and prove uniform convergence for all Levy triplets 
with this fixed value of o". To this end we show uniform convergence in 
Lemma [2j In order to control the error when going over to smaller domain 
of integration in (6.22) the term F~^{5{y + 6Y){U{u — v))fu{u)gu{v)/{uv) 



needs to be bounded uniformly. The inverse Fourier transform is bounded 
by the L^-norm of 5"^. The functions fu and gu are uniformly bounded 
since ||J^//||oo < R- The crucial step is the limit (6.23), where the refined 



dirac sequence argument is applied. As stated in (8.13), a translation before 



the Fourier transform is equal to a multiplication by a complex unit after 
the Fourier transform. Since \6\ < T{a'^^-^ + R) the complex unit tends 
uniformly to one. wu and wu do not depend on the Levy triplet. Since 
there is h with \F^{u)\ < Rh{u) and h{u) — t- as |n| — t- oo the factor 
{—u + i/U)/ exp{TFfi{Uu)) converges to —u uniformly over TisiR, CTmsL^)- 



This leads to uniform convergence in (6.23) and thus in (6.24). To see 



that the covariance without conjugation converges in (6.25) uniformly to 
zero we observe that F~^{5{y + 6)'^){u) — )• for |n| — t- oo uniformly since 
the translation by 9 corresponds to a multiplication by a complex unit by 



equation (8.13). We immediately obtain the uniform convergence in each 



component of the covariance matrix in Lemma [3} 
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8.3 Uniform convergence of the remainder term 



To begin with, we observe that Proposition [T] holds with the same constant 
for all Levy triplets in HsiR, o"max)- This follows from the bound (6.7) in the 
proof and from the fact that X{u) = e^^''^6{x)dW{x) does not depend 
on the Levy triplet. 

Lemma [5] and Lemma [6] hold uniformly, meaning that for all r/ > 



sup 

TgHs{R,(to) 

sup 



as e 



eexp(TcT2C/(e)V2) 



> 77 



0, 



e[/(e)3/2 



0. For Lemma [6] this follows from the uniform convergence of the 
bound in the proof. For Lemma [5] this can be seen by the corresponding 
uniform statements along the lines of the proof up to the bound (6.29). Then 
we bound (6.29) in two different ways depending on whether £ [0, i7g/2] 
or (T^ G ((7q/2, (Tq]. In the latter case we can proceed as in the proof for 
a > 0. Since o"^ is bounded from below by ctq/I > the convergence is 
uniform. For cr^ G [0,cjg/2] we estimate 



C 

eexp(Tc7gC/2/2) 



1 e^(;7^ + [/4)nexp(Ta2C/V 



^^dn 

T2 exp(2r(a2/2 + 7 - A) - 2T\\Tfi\\^) 



< 



T2 exp(2r(a2/2 + 7 - A) - 2T\\Tf,\\oo) ' 
which converges uniformly to zero by ctq > and by the assumption 

eC/Vlog(?7) exp(rCTgC/V2) ^ 

of Lemma[5| The maximum of the bounds is a bound that holds for all Levy 
triplets in the class. This shows the uniform version of Lemma [5} 
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